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We performed scanning tunneling spectroscopic experiments on hole-doped NdBa2Cu307_i . The 
dl /dV curves obtained at 4.2 K are asymmetric with clear peak-dip and hump structures. Energy 
derivatives of these curves show peaks at energies beyond the dip features. Highly precise full 
potential bandstructure calculations confirm a featureless electronic density of states in that energy 
region. Our results indicate that tunneling electrons couple to a collective mode in the Cu02 plane. 

PACS numbers: 71.20.-b, 74.20.Mn, 74.25.Jb, 74.50-(-r, 74.72. Bk 



The identification of the microscopic mechanism for 
pair formation in the unconventional superconductors 
is still a challenge. In conventional superconductors, 
apart from the isotope effect, tunneling experiments pro- 
vided the most direct evidence of electron-phonon in- 
teraction mediating the Cooper-pair formation [J, 
In case of high- Tc superconductors (HTSC), inelastic 
neutron scattering (INS) experiments and high resolu- 
tion angle-resolved photo emission (ARPES) experiments 
were useful to find collective modes in the range of 30-70 
meV 0, 0: S S 0, H, • Collective spin excitations have 
been suggested as the bosonic mode which interacts with 
electrons to produce a resonance peak at a wavevector 
(tt, tt) in INS spectra and a kink along the (tTjO) direction 
of the Brillouin zone in ARPES spectra [1, [1, i, 0, 
On the other hand, several ARPES results also show ev- 
idence of phonons as the collective modes 0, H, • Thus 
the assignment of the boson-mediating Cooper pair for- 
mation in HTSC is far from being settled. 

Due to the microscopic electronic inhomogeneity in 
these superconductors llj, tunneling experiments using 
scanning tunneling microscopy (STM) are expected to 
play an important role in resolving such issues being 
related to pairing. While tunneling, electrons can ex- 
cite a collective mode of a certain energy if the energy 
of the electrons is equal to that of the collective mode. 
Due to the inelastic interaction of the electrons with the 
mode, a new scattering channel is induced which results 
in a step-like feature in the d//dy curves [l^]- Ideally, 
these step-like features can be observed more easily in 
the form of a delta function in the energy derivative of 
the d//dy curves. In experiments, the peak heights and 
widths are dependent on the energy resolution, temper- 
ature and the modulation voltage For the super- 
conducting state of d-wave superconductors, Balatsky et 
al. predicted that features due to inelastic scattering of 
tunneling electrons would be observed in the d^//dV^ 
curves as satellite peaks at A -I- $7 in addition to a very 
weak peak at fi, where A is the energy gap and is 



the excitation energy [lj|. Indeed, Lee et al. have ob- 
served peaks in the second derivatives, which carry the 
signature of a bosonic mode in the hole-doped cuprate 
Bi2Sr2CaCu208-5 (Bi-2212) [li|. The tunneling spec- 
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tra of the cuprates hke YBa2Cu307_5 (Y-123), which 
contain an additional Cu-0 layer in the unit cell in the 
form of quasi-one-dimensional chains, often exhibited fea- 
tures which were not observed for Bi-2212 [l^[lB|- These 
additional features in the spectra impeded the general 
understanding of the electronic structure of the super- 
conducting state of cuprates and also masked the rele- 
vant features due to bosonic mode from the tunneling 
spectra of these cuprates. Very recently, Nietsemski et 
al. observed a bosonic mode from the tunneling spectra 
on Pro.88LaCeo.i2Cu04 (PLCCO), which is an electron- 
doped superconductor [l7[. In this work, we present 
evidence of a bosonic mode, which appears to be orig- 
inated at the Cu02 plane, from scanning tunneling spec- 
troscopy (STS) data obtained on a twinned single crystal 
of NdBa2Cu307_5 (Nd-123). Nd-123 belongs to the Y- 
123 family. They have isostructural unit cells. The ionic 
size of Nd^+ ions, which are paramagnetic at 4.2 K p^ . 
is bigger than that of Y^+ . 

The single crystal used for the measurements was 
grown from a BaO / CuO flux. The details of the growth 
process are according to those mentioned in Ref. (18] . 
After the growth, the crystal was properly oxidized to 
achieve a Tc onset of 93.5 K. The transition width of 
ATc ~ 3.5 K was determined from ac-susceptibility mea- 
surements. For the STS experiments, the as-grown crys- 
tal surface was cleaned with absolute ethanol and dried 
by pure helium gas. No special surface preparation was 
performed [4^ . The STS measurements were carried out 
at 4.2 K in He gas atmosphere using a home-made STM. 
A mechanically cut Pt-Ir tip was used for the measure- 
ments. The tunneling parameters for the measurements 
were set as -0.1 V and 0.2 nA giving a tunneling resis- 
tance of 500 Mil. The bias voltage was applied to the 
sample so that a negative (positive) voltage refers to a 
filled (empty) sample energy state. Spectroscopic data 
were obtained as [d//dy](V) curves using the standard 
lock-in technique where a modulation voltage of 2 mV 
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FIG. 1: (Color online) dJ/dV curves (shifted along y axis for 
clarity) obtained at 4.2 K along a 20 nm line. The peaks and 
humps are indicated by dotted vertical lines. 



(rms) was used. Data were taken at many different lo- 
cations on flat areas of the sample surface. The spectra 
show that the surface is electronically inhomogeneous. 
However, at some locations of the sample, we obtained 
curves which are nearly homogeneous in the length scale 
of 20-25 nm. Here we discuss only those spectra which 
contain coherence peaks (P), dips (D) and humps (H). 
These features are hallmarks of the superconducting state 
in cuprates as was independently observed in STS and 
ARPES experiments on Bi-2212 [H H^, HU, [H . In Fig. 
[U a set of d//dy spectra obtained along a line cut of 20 
nm is shown. These spectra were observed while taking 
a dl/dV curve at every 0.5 nm of a 32 x 32 nm^ area. 
It is remarkable that the PDH features are observed on 
the as-grown sample surface. On the as-prepared surface 
of NdBCO, Ting et al. found the CuO chain layer as 
the surface terminating layer [2^. The spectra shown in 
Fig. 1 do not have any similarity to the ones observed on 
the well characterized CuO chain layer [23] . It is possible 
that the topmost CuO chain layer in the present case is 
an insulating layer due to surface degradation. In that 
case, at the small bias voltage range of ±100 mV, the 
tunneling would take place between the metal tip and 
the Cu02 plane layer, which is situated within ~ 4.2 A 
of the CuO chain layer, thus revealing the clear PDH fea- 
tures. Apart from the PDH features, the asymmetry of 
the coherence peaks is evident from the data. The peak 
height is found to be always larger at the filled sample 
states compared to the empty ones. Similar asymmetric 
curves were also observed for Bi-2212 [ll|, [l^, HBj and 
Nd-123 [2^ earlier. The asymmetry is most likely re- 
lated to the probability of electron extraction and injec- 



tion from/to the material [23|. Beyond the peaks, there 
is an asymmetric V-shaped background on top of which 
humps are evident both at empty and filled states. The 
humps are observed at 2A and there are dips at ^ 
1.4 A. The respective features are broader for the empty 
states than for the filled ones. This becomes clear from 
the representative curve shown in Fig. [21 The average 
energy gap measured from peak to peak is 70 meV which 
leads to 2A/kBTc = 8.68. 

Recently, Ngai et al. reported subgap features in ad- 
dition to the coherence peaks in the dl /AV curves ob- 
tained on Yi_a;Caa;Ba2Cu307_5 thin films [iGj . Similar 
features were also reported earlier on Y-123 single crys- 
tals 15|. On the other hand, subgap features were not 
observed in the ARPES data on cleaved single crystal 
and STS data on thin film of Y-123 [H, H. However, 
the PDH features are very clearly present in those data. 
The subgap features were also not observed by Wu et al. 
on Nd-123 [26] . In cuprate systems with an additional 
copper-oxygen chain layer, the subgap features were sug- 
gested to be due to multiband superconductivity [16*1 . In 
our data, we do not find any indication of such an influ- 
ence of the metallic chain layer on the electronic structure 
representing superconducting state (see Fig. [T]). Inter- 
estingly, Derro et al. observed strong subgap resonances 
in the energy scale of ~ 5 meV while tunneling to the 
chain layer of Y-123 single crystals |2^]. These subgap 
features were suggested to be possibly due to the oxy- 
gen vacancies in the chains. It is most likely that the 
metallic CuO chain layers acquire superconductivity due 
to the proximity effect [l^, . Within this scenario, the 
spectra with subgap features possibly reflect the LDOS 
of CuO chains and those with clear PDH features reflect 
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FIG. 2: (Color online) Typical Al /AV curve (triangles) con- 
taining peak-dip-hump structures together with the energy 
derivative (circles). The smooth dark line is a guide to the 
eyes. The arrows show the inelastic features at both empty 
and filled states. The boson mode energy (fl) at empty state 
is shown. 
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the LDOS of Cu02 plane layers as these features are also 
observed on chain-less Bi-2212. If the superconductiv- 
ity in 123 systems would be due to multiband effects, 
the signature of this should also be seen in the electronic 
spectra of Nd-123. Although from the present data it 
is not clear if the subgap states would be observed in 
the LDOS of Cu02 plane that exists in the bulk, how- 
ever, it gives another indication that the PDH features 
are the generic features in the Cu02 plane-derived DOS 
of cuprate superconductors. 

In order to gain further information from the spectra, 
we focus on the part of the curves beyond the coherence 
peaks. By numerically differentiating the curves we ob- 
serve a peak in the empty state at an energy Eres as 
shown in Fig. [21 The above mentioned peaks in the 
d^I/dV'^ curves correspond to very weak step-like fea- 
tures in the dl/dV curves. The observed fine structure 
is consistent with the predicted inelastic electron tunnel 
spectroscopic (lETS) features in d-wave cuprates [isllsH. 
It also has considerable similarityto the lETS features 
observed for hole-doped Bi-2212 [14| and electron-doped 
PLCCO superconductors [17]. In order to check if the 
fine structure in the STS spectra is due to the spe- 
cific bandstructure of the superconductor, full potential 
density functional calculations were performed. We ap- 
plied the full potential local-orbital minimal basis scheme 
FPLO [13 (version 5.00-19) within the local spin den- 
sity approximation (LSDA), where the exchange cor- 
relation potential of Perdew-Wang [SS^ was chosen for 
the calculations. The strong Coulomb repulsion in the 
localized Nd 4/ shell has been modeled in a mean- 
field way within the LSDA-I-J7 approximation. A typ- 
ical value oi U4f = 8 eV has been applied throughout 
the calculations. Our results are basically independent 
of the choice of U within physically reasonable limits. 
As basis set, Nd 4/5s5p/6s6p5c?, Ba 5s5p/6s6p5d, Cu 
3s3p/4s4p3ci -I- M and O 2s2p3d + 3s states (notation: 
semicore/valence-|-polarization states) have been chosen. 
The extent of the valence states is optimized with respect 
to the total energy [IJ]. The inclusion of the semicore 
states is necessary to account for their nonnegligible over- 
lap. The polarization states provide a more complete ba- 
sis set to insure highly accurate bandstructure and DOS 
results. A very fine fc-mesh of 16200 points in the Bril- 
louin zone (2560 in the irreducible wedge) was used to 
resolve the singularities in the DOS near the Fermi level 
accurately. Convergency with respect to basis set and 
fc-mesh was carefully checked. The structural data of 
Ref. [11] were used. 

We find a total valence band width of about 9 eV, 
typical for the quasi-two-dimensional cuprates. Three 
bands cross the Fermi level, two of them originating from 
the Cu02 layer and one originating from the CuOa chains 
along the y direction. In Fig. [31 we show the calculated 
total DOS of Nd-123 in the vicinity of the Fermi level. 
In this energy region, only the states originating from 
the Cu02 planes and the CuOa chains exhibit significant 
contributions. Whereas the Cu02-planes-related DOS is 
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FIG. 3: (Color online) Calculated total electronic density of 
states (shaded gray) in the vicinity of the Fermi level (zero 
energy) for NdBCO. The sharp features at about 85 meV and 
110 meV binding energy can be assigned to van-Hove singu- 
larities originating from the Cu-O chains ([red] dot-dashed 
line). The contribution from the plane ([green] dashed line) 
is rather flat in the shown energy window. 



almost constant close to the Fermi level, two distinct van 
Hove singularities related to the CuOs chain show up at 
about 85 and 110 meV. Otherwise, no significant features 
in the DOS near the Fermi level are obtained. The slight 
underdoping of the present sample can be modeled by 
a small rigid shift of the Fermi level towards negative 
energies. For an oxygen deficiency S = 0.01, this shift is 
estimated to be -18 meV. There is no sharp structure 
present in the DOS of the Cu02 plane at energies close to 
the region where peaks are observed in the experimental 
d?I /dV'^ spectra. Thus, these features cannot be related 
to the DOS. The calculations also show that the bilayer 
splitting is not responsible for the PDH features. 

In accordance with Eliashberg's classical theory for 
strong coupling superconductors, where the electron- 
boson interaction was proposed to be observed at an en- 
ergy oi E — IS. + VL [1^ , the VL values were determined for 
both empty and filled state. The empty and filled state VI 
are experimentally determined by ^{r) = Eresij") — '^(f) 
from the entire set of curves obtained at different loca- 
tions (r) on the sample. The mean values of f2(r) are 
found to be 22.9 ± 1.8 mcV and 23.7 ± 1.6 meV for 
empty and filled states, respectively. Thus, the peaks at 
Eres in the d'^I/dV'^ curves are found to be almost sym- 
metric. This strongly supports the assumption that the 
observed features are signatures of inelastic tunneling. 
They result from an additional tunneling channel intro- 
duced due to the excitation of a collective mode at ~ 
23 meV. The mode energy is surprisingly low compared 
with the ones observed for Y-123 for similar dopings [3]. 
A priori, it is not clear if the mode is excited in the tun- 
neling barrier or in the Cu02 plane of the sample. Thus, 
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FIG. 4: (Color online) Mode energy (^(r)) vs energy gap 
(2A(r)) obtained from the dl/dV curves like in Fig. 2. The 
dashed line emphasizes the tentative correlation. The uncer- 
tainties in the evaluation of A and are ±(0.5 • • ■ 1) meV and 
±(1 • ■ ■ 1.5) meV, respectively. 



it is instructive to plot ri(r) against 2A(r). As shown in 
Fig. m il and A are intimately related. A similar, but 
comparatively strorigerj correlation was also found for Bi- 
2212 and PLCCO UdUJl. Moreover, this mode energy 
for Nd-123 was observed for only those spectra which ex- 
hibit clear PDH features. Thus, the collective mode ap- 
pears to be intrinsic to the superconducting Cu02 plane. 
To determine the nature of the observed collective mode, 
we compare the mode energy with the known phonon en- 
ergies of Nd-123. The former is rather low compared to 
the in-phase or out-of-phase 0(2)-0(3) vibrations which 
have energies of more than 35 meV as observed from Ra- 
man spectra [13, [H, [s^. In case of Bi-2212, Pilgram 
et al. suggested that the peaks in the second derivative 
are due to the excitation of an apical oxygen vibration 
mode [4C|. However, the corresponding energy for Nd- 
123 as detected from Raman spectra is ^ 68 meV [37| . 
Also, this is not intrinsic to the Cu02 plane. The only 



low-energy modes for Nd-123 are due to the metallic Ba 
(14 meV) and Cu (19 meV) ion vibration modes. The 
latter one seems to be the only possible mode that could 
be excited during tunneling in the present case. However, 
it is not clear why in that case the same mode energy is 
not observed in case of Bi-2212 or PLCCO [lllI3|. Thus, 
although the mode does not appear to be a phonon, we 
can not definitely rule out the possibility of a low energy 
excitation of metal ions. 

For an underdoped Y-123 sample, incommensurate 
spin fluctuations were observed at 24 meV in INS ex- 
periments [4l|. Niestemski et al. have suggested the 10.5 
meV mode as observed for PLCCO to be a magnetic res- 
onance mode [l3|. At this point, it remains unclear what 
the origin of the observed mode for this rare-earth-based 
cuprate is. Since the features observed in the spectra and 
the correlation between A and O are qualitatively simi- 
lar to those observed for Bi-2212, the identification of the 
nature of this low energy mode in the rare-earth-based 
cuprate would be important. 

In summary, our STS data on as-grown slightly under- 
doped Nd-123 single crystal do not show any evidence 
of multiband superconductivity. However, it would be 
important to probe the Cu02 plane layer in the bulk 
in order to completely rule out the role of CuO chain 
band towards superconductivity in these systems. The 
peak-dip-hump features in the spectra indicate that these 
are most likely the generic features of the Cu02 plane- 
derived LDOS in cuprates. Furthermore, a bosonic mode 
has been detected with a characteristic mode energy of 
about 23 meV. The mode energy is rather low compared 
to the other bosonic modes observed for the cuprates. 
The observed mode energy has certain correlation with 
the energy gap and thus appears to be intrinsic to Cu02 
plane. However, the correlation is weak compared to the 
ones reported for other cuprates. Thus, the origin of the 
low energy mode due to a tunneling pathway, which is 
extrinsic to the Cu02 plane layer, can not be completely 
ruled out. 
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